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We study in detail the d = 6 operators for proton decay in the two possible matter 
unification scenarios based on SU(5) gauge symmetry. We investigate the way to distinguish 
between these two scenarios. The dependence of the branching ratios for the two body 
decays on the fermion mixing is presented in both cases. We point out the possibility to 
make a clear test of flipped SU(5) through the decay channel p — > tt + D, and the ratio 
T (p^K°e+)/T(p^ir°e+). 

I. INTRODUCTION 



Proton decay d] is the most dramatic prediction of grand unified theories, where quarks and 
leptons are at least partially unified. Its signatures have been extensively studied in various the- 

for many years. Recently, in the context of minimal 
supersymmetric SU(5), the predictions coming from both d = 5 and d = 6 operators have been 
studied in order to understand if this model is ruled out Several solutions have been 

forwarded 13|, Il4t lla ] to address this issue. This has renewed the interests of many groups to the 
important question of the proton stability (for a review see ^(j]). Similar study, in the context of 
flipped SU(5), has also been made jl^] concluding that the flipped model is out of trouble. 
There are several contributions to the decay of the proton. The d = 4 and d = 5 are the most 
important in supersymmetric scenarios. In a theory where matter-parity is conserved the d = 4 
are forbidden, while the d = 5 operators can always be suppressed by choosing a particular Higgs 

are the d = 6, which we study here in 

detail. 

An extensive study of d = 6 operators in the most general way in the context of SU(5) and SO(W) 



has been preformed in reference 



2l| . There, it has been pointed out that it is possible to make a 



clear test of any grand unified theory with symmetric Yukawa couplings through the decay chan- 
nels into antineutrinos. However, the particular case of flipped SU (5) has not been studied taking 



21124)- With 



into account the general dependence on fermion mixing (for early analyses see |22l . 
this work we seek to remedy that. Namely we investigate all d = 6 proton decay operators in 
two different GUT models based on SU(5). We then confront the signatures of the two unifying 
schemes pointing out the way to distinguish between them. We also point out the way to test 
flipped SU(5). 

The manuscript is organized as follows. In Section [H] we briefly review the key properties of both 
SU(5) and flipped SU(5) unified theory. Section ITTT1 is devoted to the general discussion of d = 6 
operators in both scenarios. In Section llVl we specify all the branching ratios for the independent 
channels for proton decay. That section contains the main results of our work. Finally, we con- 
clude in the last section. The Appendix contains useful decay rate formulas used throughout the 
manuscript. 



II. MATTER UNIFICATION BASED ON SU(5) 

The smallest special unitary group that contains the Standard Model (SM) gauge group is 
SU(5). The SU(5) grand unified theory 2^, [2^ is an anomaly free theory, where we have partial 
matter unification for each family in three representation: 10, 5 and 1. The singlet is identify with 
the right handed neutrino. In the SM language we have: 10 = (3, 2, 1/3) © (3, 1, -4/3) © (1, 1, 2) = 
(Q,u c ,e c ), 5 = (3,1,2/3) © (1,2,-1) = (d c ,L), and 1 = (1,1,0) = v° ', where Q = (u,d) and 
L = (v, e). The off-diagonal part of the gauge fields residing in the 24 of SU(5) is composed of 
bosons (X, Y) = (3, 2, 5/3) and their conjugates, which mediate proton decay. X and Y fields have 
electric charge 4/3 and 1/3, respectively. 

The electric charge is a generator of conventional SU(5). However, it is possible to embed the 
electric charge in such a manner that it is a linear combination of the generators operating in both 



SU(5) and an extra U(l 
done in a flipped SU (5) 



and still reproduce the SM charge assignment. This is exactly what is 
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27, 



28, 29]. The matter now unifies in a different manner, which can 
be obtained from the SU(5) assignment by a flip: d c *-* u , e c <-> v , u <-> d and v «-> e. In the 
case of flipped SU(5) the gauge bosons responsible for proton decay are: (X',Y r ) = (3,2, —1/3). 
The electric charge of Y' is —2/3, while X' has the same charge as Y. Since the gauge sector and 
the matter unification differ from SU(5) case, the proton decay predictions are also different ^22]- 
Flipped SU(5) is well motivated from string theory scenarios, since we do not need large repre- 
sentations to achieve the GUT symmetry breaking j^. Another nice feature of flipped SU(5) is 
that the dangerous d = 5 operators are suppressed due to an extremely economical missing partner 
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mechanism. This allows us to concentrate our attention to the gauge d = 6 contributions. 
We next analyze the possibility to test two realistic grand unified theories: the SU(5) and flipped 
SU(5) theory. We make an analysis of the operators in each theory, and study the physical pa- 
rameters entering in the predictions for proton decay. We do not commit to any particular model 
for fermion masses, in order to be sure that we can test the grand unification idea. 



III. D=6 OPERATORS 



In the Georgi-Glashow SU (5 
to the decay of the proton are 



matter unification case, the gauge d = 6 operators contributing 



(la) 
(lb) 



®SU(5) — ^1 ^' l jk U< ^a. 7^ Qj aa 7a* Qkfibi 

°SU(5) = k l e Vk U ?a 7^ Qjaa dg 7/1 L ph . 



On the other hand, flipped SU (5) matter unification yields: 

SU(5)' = k 2 e ijk e «/3 dg 7^ Qjp a ug ^ L ab , (2a) 

°SU(5)' = k 2 ei 3 k e °/3 d< ta 7^ QjPa 7jtt Qkab- (2b) 

In the above expressions k\ = g^M7 XY y an d ^2 = H^MT^, Y /y where M(x,Y) (M(X',Y')) ~ Mqut ~ 
10 16 GeV and c/5 (g' 5 ) are the masses of the superheavy gauge bosons and the couplings at the GUT 
scale in SU(5) (flipped SU(5)) case, i, j and k are the color indices, a and b are the family indices, 
and a, j3 = 1, 2. 

In these theories the diagonalization of the Yukawa matrices is given by the following bi-unitary 
transformations: 



U T C Y V U = 
D T C Y D D = Y^\ 



E T C Y E E 



v diag 
E ■ 



(3) 
(4) 
(5) 



Using the operators listed in Eqs. (|T|). the effective operators for each decay channel in the SU(5) 



2l|: 



case upon Fierz transformation take the following form in the physical basis 

0(e% , dp)su(5) = c(e%,dp)su(5) djk up ^ uj e% 7 M d kj3 , 

0(e a ,d^) SU (5) = c(e a ,d^) su{5) e ijk uf 7^ Uj d^ 7^ e a , 

0(v h d a ,dp) SU ( 5 ) = c(vn,d a ,djj) S u(5) e ijk u ? 7 M d ja d%p 7 M v h 

0(uf,d a ,d^) su ^ = c(vf,d a ,d$) su( s) e ijk dfy 7^ uj uf 7 M d ka 



(6a) 
(6b) 
(6c) 
(6d) 



4 



where 



c(ea,d/3)su(5) 
c(e a ,dp)su{S) 

c(p t ,d a ,dp)su(S) 



k 



0. 



In the case of flipped SU(5) (see Eqs. @) the effective operators are 



0{e%,dp) SU ( 5 y 



0(e a ,d c 



■0 )SU(5)' 



0(vi,d a ,dp)su(b) 
0(v>i ,d a ,d /3 )su(5) 



c{ea,df3)su(5Y tijk uf 7 M Uj eg 7/Lt 4^, 
c(e a ,^) 5[/{5) / e iifc 7^ uj d%p 7^ e a , 
c iyhd a ,dp)su(^y €ijk uf 7^ 
c(vf,d a ,dp) SU ( 5 y e ijk dfp 7 M u 3 - i/p 7/J d fccr , 



where 



c ( e aidf3)sU(5)' 

c(e a ,dp)su(5)' 
c(vi,d a ,dp)su(5)' 
c(vi ,d a ,dp)su(5)' 



0, 



> lo 



fefVf "(^iVaD^VsV^jv)", « = 1 or p = 1 

(v 4 vt ?\u EN v 2 r + 



^"2 



o 



1 or 



(7a) 
(7b) 
(7c) 
(7d) 



(8a) 
(8b) 
(8c) 
(8d) 



(9a) 
(9b) 
(9c) 
1. (9d) 



We use the subscripts SU(5) and SU(5)' to distinguish the two scenarios. The mixing matrices 
Vx = U j c U, V 2 = E j c D, V 3 = D j c E, V A = D f c D, V UD = U^D, Ven = E^N and U EN = E j c N c . 
The quark mixing is given by Vjjd = U^D = K\VckmK2, where K\ and K2 are diagonal matrices 
containing three and two phases, respectively. The leptonic mixing Ven = K^V^K^ in case of 
Dirac neutrino, or Ven = K3V1 in the Majorana case. Vj° and V t M are the leptonic mixing 
matrices at low energy in the Dirac and Majorana case, respectively. 

Notice that in general to predict the lifetime of the proton in SU(5), due to the presence of d = 6 
operators, we have to know k\, V 1 , V2, V3, while in flipped SU(5) we have to know &2, V\ , V3, 
V4 and Uen- In addition we have to know three diagonal matrices containing CP violating phases, 
K±, Ki and K3, in the case that the neutrino is Majorana. In the Dirac case there is an extra 
matrix with two more phases. 

From the above equations, we see that there are no decays into v in SU(5), and in flipped SU(5) 
into e , since these are singlets in the corresponding scenarios. 
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IV. FLIPPED SU(5) VERSUS SU(5) 

There are only seven independent relations for all coefficients of the gauge d=6 operators con- 
tributing to nucleon decay 2l|. Therefore, if we want to test a grand unified theory, the number 
of physical quantities entering in the proton decay amplitude must be less than that. This is im- 
portant to know in order to see if it is possible to test a GUT scenario. 

Since we cannot distinguish between the neutrino flavors in the proton decay experiments, in order 
to compute the branching ratios into antineutrinos we have to sum over all of them. Using the 
expressions in the Appendix, and the following relations: 



Y,c{vud a ,d% u{b) c{v u d^d c & ) su{b) = kf ( Vf* Vu D ) la (ViVud) 17 ^ 135 , (10a) 
1=1 

3 

Y,c^i,d a ,d^r su{5Y c^,d^df) su{5 y = ki(v:f a v^, (iob) 
i=i 

we can write down the ratios between the lifetimes in both theories for the decays into antineutrinos. 
They are: 



(11a) 
(lib) 
(11c) 



r(p - 




r(p - 




r{p - 


- 7T+u) su ( 5 y 


r{p - 




r(n - 


-» K°v) su W 


r(n - 





k\ 


AIKVx^Vckm) 11 


2 + Al 


{VxKxVckm) 12 


2 




y, 21 \ 


2 + A\\vp 


2 + A x A 2 {{VlYWP + (Vf) 


12^21) 



kf KViKxVckm) 
h o 



11 12 



2 \V, 



1 



kf 


AlKVtKxVcKAi) 11 


2 + Al 


(V^Vckm) 12 


2 




n 21 l 


2 + Al\vP 


2 + A 3 A 2 ((U 4 *) 21 U 4 12 + (V 4 *) 


12 F 4 21 ) 



where 



2m n 

M = - — -D, 12a 

Stub 

777 

A 2 = 1 + -JL(D + 3F), (12b) 
6m b 

A 3 = 1 + ^(D-3F). (12c) 
6m b 



The same procedure can be done for the decays into charged leptons: 

2 



r(p -> vr°e+) 5C/ ( 5 ) 
r(p -> K°e+)^( 5 ) 



+ 



k-2 



{V 4 K*V£ KM K* ) 11 (V! K x V CK mK 2 V\ V 3 ) ^ 



-(13) 



+ 



Vl x v£ 2 + iV 1 K 1 V CK MK 2 ) 12 iV 2 K^ KM K*f 1 



k'2 



{V 4 K*vi KM K*)^iV 1 K 1 VcKMK 2 vlv 3 )^ 



-(14) 



Eqs. (|llj). (|13|). and (fTl]) are the most general equations that we could write in the two scenarios 
and will help in future to distinguish between them if proton decay is found. In other words, for 
a given model of fermion masses, using the above equations we could see the difference in the 
predictions for proton decay. Unfortunately, as one can appreciate, the branching ratios depend on 
too many unknown factors, including the new CP violating phases. (These, in principle, could be 
defined in a particular model for CP violation.) Therefore it is impossible to test those scenarios in 
general through the decay of the proton unless we known the flavor structure of the SM fermions. 
Since we cannot make clear predictions in the most general case, let us consider special cases in 
these two matter unification scenarios based on SU(5) and compare them. 



A. SU(5) with Y v = Yff 

In SUib), if Yjj = Yjj , we have Uc = UK U , where K u is a diagonal matrix containing three CP 
violating phases. Therefore we get: 

3 

Y,<n,d a ^r su{ , )C { n ,d^) su ^ = kKvs KM ) ia iK* 2 r a ivcKMr'K^5^. (15) 
i=i 

In this case, as has been shown j^ . the clean channels, i.e., the channels that we have to look to 
test this scenario, are: 



T(p^K+D) = k{[A\\Vl, l KM \ 2 + Al\Vh 2 KM \ 2 ]Ci, 
Tip^K+v) = kf |^ M | 2 C 2 , 



(16a) 
(16b) 



where 



C 2 



[mz 



m 



2 \2 

k) 



A? I |2 

-A L \a\ 



E-Al\a\ 2 il + D + F)< 



(17a) 
(17b) 
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Notice that we have two expressions for k\ , which are independent of the unknown mixing matrices 
and the CP violating phases. Therefore it is possible to test SU(5) grand unified theory with 



symmetric up Yukawa matrices through these two channels |2l( . Notice that these results are valid 
for any unified model based on SU(5) with Yjj = . For example, this includes the case of minimal 
SUSY SU(5) with two extra Higgses in the fundamental and antifundamental representations. The 
case of modified missing doublet SUSY SU (5) model 0, 0| is also included in our analysis. 



B. Renormalizable flipped SU(5) 

In renormalizable flipped SU(5) we have Yd = , so Dc = DKj, where is a diagonal 
matrix containing three CP violating phases. In this case the coefficients entering in the proton 
decay predictions are: 



(18a) 



i=i 



c(e a , (ijg ) | — /c 2 



V, 



1/3 



CRM 



(ViV UD V}v 3 ) 



la 



A'2 



V, 



CRM 



(18b) 



Using these equations we get the following relations: 



r(p-vr°4 

T{p -> K°e+ 

r(p - vrO e +; 



- r(p 

2 yF 



7T+P) W, 



ii 

CRM] 



(u c Ey 



2^3 rCRM 



c 2 \v n 

z VCRM 



2 • 



(19a) 
(19b) 
(19c) 



where: 



(m; 



m 



2 \2 
R) 



-A\\a\ 



1 + HhL( D _ F ) 

m B 



(20) 



Notice that in this case, T(p -> #+z?) = 0, and T(n -> # Q z?) = 0. In Eq. (FTT^) we assume 
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(U^E) 1 * ± 0. 

We can say that the renormalizable flipped SU(5) can be verified by looking at the channel 
p — > 7r + P, and using the correlation stemming from Eq. (|19cl) . This is a nontrivial result and 
can help us to test this scenario, if proton decay is found in the next generation of experiments. It 
is one of the main results of this work. If this channel is measured, we can know the predictions 
for decays into charged leptons using Eq. (|19b|) for a given model for fermion masses. Therefore it 
is possible to differentiate between different fermion mass models. 

Note the difference between Eqs. ()16b|) and (|19a|) : there appears a suppression factor for the chan- 
nel p — > tt + V in the case of SU(5). 

Since the nucleon decays into K mesons are absent in the case of flipped SU(5), that is an inde- 
pendent way to distinguish this model from SU(5), where these channels are always present. 

V. CONCLUSIONS 

We have investigated in model independent way the predictions coming from the gauge d = 6 
operators in the two possible matter unification scenarios based on SU (5) gauge symmetry. We 
write down the most general ratios between the lifetimes in SU(5) and flipped SU(5) theory for 
each channel, providing the way to distinguish between them. We find that in general it is very 
difficult to test flipped SU(5). However, in the case of renormalizable flipped SU(5) model, the 
decay channel p — > tt + 17, which is a clean channel, and the ratio r(p — > K°e^)/T(p — > vr°e+) could 
be used to test this theory. If the decay of the proton is found in future, our results will be useful 
to analyze the predictions in these theories. 
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APPENDIX 

Using the chiral Lagrangian techniques (see reference [3]), the decay rate of the different 
channels due to the presence of the gauge d = 6 operators are given by: 



r(p -» i^ + z 



r(p -» vr+p) 

r(p -» T/eJ) 

r(p - i^°e+) 
r(p-vr°e+) 

r(n -> if°77) 



r — > (n — ► 77^) 
r(n -» vr~e+) 



(m 2 -mf^) 2 



8vrm3/2 
3 



A| |a| 2 (1 + D + F) 2 | C K d > d 



8=1 



(m 2 - m 2 ) 2 
48vr/ 2 m3 



cfa, d, s C ) + [l + pL( D + 3F)]c(u t ,s, d c ] 
3m b am b 



C\\ 2 



i=l 



»? ; 42 1 ,|2 



A£ |a| a (1 + D - 3F) 2 {|c( e/3 , d c )\* + |c(e£, d)D 



-C 



(m 2 - m 2 K ) 2 



A\ \a\ 2 [1 + —?-(D - F)] 2 {|c( e/3 , + \c^, a )\'} 



m 



CM 2 



c 



87r/ 2 m 3 



m B 



m r , 



^A 2 L \a\ 2 (1 + D + F) 2 {\c{e p , d c )\* + |c(eg, d)| z } 



2 \2 



16tt/, 

(m 2 - m^ ) 
87rm3/2 

3 

8=1 

m- 



(A.l) 
(A.2) 
(A.3) 

(A.4) 
(A.5) 



2 U,|2 



-Af. la 



c(ui,d, s c )[l + ^(D- 3F)] - c(v i: s, d c )[l + ^(D + 3F)} 



3m B 



3m b 



(A.6) 



T(n^7T°u) = -^A 2 L \a\ 2 (l + D + F) 2 Y,\c(^,d,d 



CM 2 



16VT/2 

(m 2 - m 2 ) 2 
48vrm3/2 



i=l 



A| H 2 (1 + U - 3F) 2 ^ |c(fj, d, d 



m 



n A\ |a| 2 (1 + D + F) 2 {\c(ep, d c )\* + \c(e%,d)\'} 



8=1 

'C\ " 



8vr/ 2 



(A.7) 

(A.8) 
(A.9) 



In the above equations mg is an average Baryon mass satisfying ~ my, ~ roA, -D, i 7 and a are 
the parameters of the chiral lagrangian, and all other notation follows [32(. Here all coefficients 
of four-fermion operators are evaluated at Mz scale. Al takes into account renormalization from 
Mz to 1 GeV. Vi = u e , v T and ep = e, //. 
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